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Abstract The core structure of galaxy clusters is fundamentally important. Even 
though self-gravitating systems have no stable equilibrium state due to their nega- 
tive heat capacity, numerical simulations find density profiles which are universal in 
the sense that they are fairly flat within a scale radius and gradually steepen farther 
outward, asymptotically approaching a logarithmic slope of « —3 near the virial ra- 
dius. We argue that the reason for the formation of this profile is not satisfactorily 
understood. The ratio between the virial radius and the scale radius, the so-called con- 
centration, is found in simulations to be closely related to the mass and the redshift 
and low for cluster-sized haloes, but observed to be substantially higher at least in a 
subset of observed clusters. Haloes formed from cold dark matter should furthermore 
be richly substructured. We review theoretical and observational aspects of cluster 
cores here, discuss modifications by baryonic physics and observables that can provide 
better insight into the internal structure of clusters. 



1 Fundamental importance of cluster cores 

1.1 Difficulties with self-gravitating systems 

The core structure of galaxy clusters is of fundamental importance not only for as- 
trophysics, but for physics as a whole. The essential reason for this statement is that 
the properties of the yet hypothetical dark-matter particles, and the dominance of not 
electromagnetically interacting dark matter itself, can be tested in the cores of galaxy 
clusters. 

Due to the virial theorem, self-gravitating systems have a negative heat capacity. 
This is easily seen as follows. In its scalar version, the virial theorem demands that the 
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total kinetic energy i?ki n and the potential energy E po t be related by 

2-E kin = -Spot • (1) 

The factors multiplying the two forms of energy are due to the fact that the kinetic 
energy is a homogeneous function of degree 2 of velocity, while the potential energy is 
a homogeneous function of degree —1 of position. The total energy is thus 

E = -Ekin + Spot = —Skin . (2) 

Self-gravitating systems can lose energy in various ways. For example, particles can 
become unbound and be ejected from the system in three-body encounters, or baryonic 
gas mixed with the dark matter can radiate energy and adiabatically contract the 
dark matter. If a self-gravitating system loses energy, no matter how, the system's 
total energy decreases. Since it is negative for a bound system, its absolute value 
increases thereby and, following the virial theorem (2), its kinetic energy must also 
increase. Cooling of the system by any kind of energy loss thus leads to heating of 
the system by contraction, which defines a negative heat capacity. This in turn implies 
that self-gravitating systems cannot reach a stable equilibrium state where a suitable 
thermodynamical potential could attain a minimum. 

Nonetheless, numerical simulations of dissipation-less dark matter under the ex- 
clusive influence of gravity unanimously show that long-lived, gravitationally bound 
systems form in a wide range of masses, the so-called dark-matter haloes. Quite inde- 
pendent of various initial conditions such as the power spectrum of the dark matter 
fluctuations, these haloes turn out to have a universal density profile p(r), characterised 
by two radii. Within the outer radius, the virial radius 7? v i r to be defined below, virial 
equilibrium is supposed to be achieved. Qualitatively speaking, the density profile is 
steep from an inner radius, the scale radius r s , to the virial radius, and flat from the 
halo centre to the scale radius. While simulations agree that the double-logarithmic 
slope 

dmr 

reaches the asymptotic value —3 near the virial radius, inner slopes between —1 and 
— 1.5 are typically measured in simulated haloes. Introducing the dimension-less radius 
x = r/r s , density profiles of the functional form 

p(x) = — - Ps ,, (4) 

' y ' x a {l + x) 3 - a v ' 

fit the numerically simulated dark-matter haloes very well, with — 1 < a < —1.5. 
For a = —1, the profile (4) was first suggested by (Navarro et al. 1997, 1996, hereafter 
abbreviated by NFW), while Moore et al. (1998, 1999) favoured steeper central profiles 
with a « —1.4. More recent, highly resolved numerical studies find that the density- 
profile slope decreases gently from the scale radius inward (Power et al. 2003; Navarro 
et al. 2004) in a way that may depend on the halo mass (Jing and Suto 2000). A density 
profile with an alternative functional form, 

/ 9(a;) = p s exp[-2n( a; 1/ ™-l)j , (5) 

was suggested by Einasto (see Einasto and Haud 1989) and found to fit density profiles 
of numerically simulated dark-matter haloes very well (Merritt et al. 2006) with n « 5. 
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Fig. 1 The NFW and Einasto density profiles (the latter with n = 5) are shown as functions 
of the dimension-less radial coordinate x = r/r B . The difference are subtle, but the Einasto 
profile's central slope gradually decreases to zero as the centre is approached. The profiles are 
scaled such that the density is unity at r = r B . 

In (5), the scale radius r s was introduced such that the profile slope is a = —2 at 
r = r s , which is also the case for the NFW profile at x = I . 
The slope of the Einasto profile, 



decreases continuously towards the centre and reaches zero there. One important dif- 
ference between the generalised NFW profile (4) and the Einasto profile (5), shown in 
Fig. 1, is the central density cusp in the former, but not in the latter. If the dark-matter 
particles decay or annihilate, the annihilation rates in both types of density profile may 
be substantially different. As two-body processes, annihilation rates are proportional 
to the squared particle density, hence cusped density profiles may predict higher rates 
than profiles with a finite-density core. More important than the central cusp, however, 
seems to be the level of substructure in cluster haloes, as we shall discuss later. 

We shall return further below to the issue of the slope a near the centre of cluster 
haloes. For the purpose of our present discussion, the exact value of a is less important, 
however. Since stable final states of self-gravitating systems cannot be expected for the 
reasons given above, the mere fact that long-lived transient states like dark-matter 
haloes do exist at all and form a universal density profile is a most remarkable finding 
on its own. 

The ratio between the two radii, the concentration parameter 



is found to be closely related to the masses and the formation times of the haloes. 
Generally, the earlier a halo forms, the larger the concentration parameter is. This 
is commonly attributed to the higher cosmic background density at the time of halo 
formation. We shall discuss the concentrations in more detail below. 



a = — 
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1.2 Hierarchical formation of dark-matter haloes 

The global properties of a dark-matter halo, i.e. its radius ii v j r and its mass M, are com- 
monly related to each other and to the parameters of the cosmic background through 
the spherical collapse model. Specifically, the model states that the mean matter den- 
sity enclosed by the virial radius should be higher than the mean cosmic density p by 
the virial overdensity Z\ v i r , thus 



3M 



= A, 



(8) 



The virial overdensity A v ^ depends on redshift and on cosmology. In terms of the 
matter-density parameter i? m rj at the present epoch or of the matter-density O m (z) as 
a function of redshift, the mean cosmic matter density is 



8ttG 



3H 2 (z) 
8nG 



(9) 



For an Einstein-de Sitter universe with critical matter density and without a cosmolog- 
ical constant, the mean overdensity is analytically found to be Zi v ir = 187T 2 ~ 178. For 
the standard cosmological model, A v (relative to the mean cosmic density) is substan- 
tially higher and changes with redshift, as shown in Fig. 2. Often, A v is instead given 
relative to the critical cosmic density and thus reduced by a factor fi m (z) compared 
to the values shown in Fig. 2. 




redshift z 

Fig. 2 The virial overdensity A v relative to the mean cosmic density is shown here as a 
function of redshift for the standard cosmological model (J? m = 0.27, Qa = 1 — &m)- 



The considerations so far show that the formation and in particular the existence 
of long-lived transient states of dark-matter haloes, at least in numerical simulations of 
collision-less dark matter, seems to be a remarkably well-determined process without 
much freedom. The collapse dynamics sets the virial overdensity A V - 1T which, together 
with the cosmic mean matter density p, relates the virial radius to the mass. The 
concentration, which itself seems to be defined by the formation redshift, then sets the 
scale radius. Given that, the density profile is defined. 
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If the dark matter is cold, the power spectrum P(k) of the initial density fluctu- 
ations from which the haloes ultimately form has a characteristic shape that is again 
fully determined by few cosmological conditions. For small k, it increases as a power 
law in k, P(k) oc k n with n < 1. For very large k, it decreases asymptotically like 
P(k) oc fc n ~ 4 . The power-law index n is set by cosmological inflation very early during 
the evolution of the universe. If inflation had lasted eternally, the power spectrum would 
be precisely scale-free, n = 1. Since inflation must have ended after a finite time, n 
should be slightly smaller. It belongs to the most exciting results of the WMAP satellite 
mission that n is indeed slightly, but significantly smaller than unity, n — 0.968 ±0.012 
(Komatsu et al. 2011). 

The decrease of the slope by —4 is also entirely defined by the inhibition of growth of 
dark-matter fluctuations inside the horizon during the radiation-dominated epoch. The 
maximum of the power spectrum is thus set by the horizon size at the epoch of matter- 
radiation equality, which is in turn defined by the cosmic matter-density parameter. It is 
essential for our purposes to note that the dark-matter fluctuation power spectrum has 
no free parameters except for its overall normalisation. Traditionally, this normalisation 
is fixed by the parameter as, generally defined by specialising the scale-dependent 
variance 

« 2 R=j Q O ° ! 0P(k)\W R (k)f (10) 

to the radius R — 8 h^ 1 Mpc. In (10), Wn(k) is a Fourier-transformed window function 
introduced to smooth fluctuations on scales smaller than R. The specific setting of this 
length scale to 8/i _1 Mpc has historical reasons: First measurements of ag based on 
the galaxy distribution found that variance of the density contrast was approximately 
unity on a scale of 8 ft -1 Mpc (cf. Davis and Peebles 1983, for an example). 

The power in the dark-matter fluctuations is proportional to the volume in fc-space 
times the power spectrum, or oc k 3 P(k), which asymptotically tends to fc" _1 for large 
k. Since n < 1, the power keeps rising with k until it gently flattens off asymptotically 
for very large k. This means that the largest power is in the smallest fluctuations, giving 
rise to the bottom-up scenario of cosmological structure formation in cold dark matter: 
Those fluctuations with the largest power form first, hence structures in cold dark 
matter start forming at the smallest scales. It is thus expected that cold dark matter 
clumps on all scales. Moreover, the largest structures will form the latest. Depending 
on the normalisation parameter ag, a maximum mass is expected which could just have 
formed in the cosmic present. 

The statistics of cosmological structures is characterised by the assumption that 
they form from an initially Gaussian random field with a low amplitude. During the 
linear phase of structure growth, this amplitude increased proportional to the linear 
growth factor D+(a) without any interaction between different Fourier modes of the 
density field. The spherical-collapse model predicts that non-linear collapse concludes 
when the linear density contrast has reached a critical value 8 C which is 



in an Einstein-de Sitter universe and not much different in other model universes. The 
abundance of cosmological objects of a certain mass can then be derived directly from 
the peak statistics in a Gaussian random field. The result is the Press-Schechter mass 
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function, 

N{Ml a) = Jl ^/^exp (- AttI ^ • (12) 

where po is the mean cosmic matter density today. The radius 7? and the mass M are 
related by requiring that a sphere with radius R filled with the mean cosmic matter 
density should contain the mass M, 

^-R 3 P = M. (13) 

Again, the details of this mass function and its possible modifications are not relevant 
for our discussion. It is important, however, that the mass function cuts off exponen- 
tially at a mass scale defined by 

D+(a)a R ^6 c . (14) 

The larger R is, and therefore the larger the masses are that are being considered, the 
lower is the variance a R and the larger D+(a) has to become for (14) to be satisfied. 
This quantifies the earlier statement that the most massive objects form late in cosmic 
history. 
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Fig. 3 Left panel: The Press-Schechter mass function is shown as a function of mass for three 
different rcdshifts. The number density of massive clusters is evolving very rapidly as a function 
of time, showing that clusters are a dynamical population of cosmic objects. Right panel: The 
non-linear mass M* is shown as a function of redshift for a standard ylCDM cosmological 
model. Even at redshift zero, the non-linear mass is just above 10 13 Mq, showing that clusters 
with their masses > 10 14 Mq are highly non-linear objects. 



A typical non-linear mass scale is commonly defined by combining Eqs. (14) and 
(10). Define a window function W R (r), for example a top-hat function Wn(r) — &(R — 
r). Its Fourier transform 

inserted into (10), returns the variance a R on the scale R. As 7? increases, progressively 
more structure is smoothed out, so a R decreases. Find the scale where ur = <5 C and 
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convert this to a mass through Eq. (13), which is called the non- linear mass M*. The 
non-linear mass for a standard /1CDM cosmological model is shown in the right panel 
of Fig. 3. 

1.3 Expected level of substructuring 

Numerical simulations at high resolution, assuming cold dark matter, find strongly 
substructured haloes. Numerous different simulation runs undertaken during recent 
years arrive at the following statements on the substructure level: 

The amount of substructure in a halo increases with the mass of the host halo. At 
fixed host-halo mass, it decreases with the concentration of the host halo and with its 
formation redshift. These trends are noticeable, but of comparable amplitude as the 
scatter in the subhalo abundance between similar host haloes. The abundance of low- 
mass substructures per unit host-halo mass reflects that of the surrounding Universe 
as a whole. The subhaloes are substantially less spatially concentrated than the dark 
matter. Different approaches to the identification of subhaloes within host haloes lead 
to similar substructure abundance functions (Gao et al. 2004; Zentner et al. 2005; 
Giocoli et al. 2010b). 

Taking into account baryons in addition to the dark matter, Dolag et al. (2009) 
find that the diffuse hot gas is efficiently removed from subhaloes when they enter their 
host haloes. Compared to simulations with purely dark matter, stripping of the gas 
may decrease the subhalo mass function. However, if radiative cooling of the baryons 
is taken into account, baryons condense towards the centre of the subhaloes and form 
compact stellar cores, thus mitigating the effect of stripping. In all variants of the 
simulations including baryons, only about 1 % of the subhaloes within the virial radii 
of their host haloes preserve a gravitationally bound hot gaseous atmosphere. 

Based on the Millennium-II simulation data, Boylan-Kolchin et al. (2009) proposed 
a fitting formula for the abundance function of subhaloes, 



where /i — (M su b/M2oo) is the ratio between the masses of subhaloes and the virial 
mass of their host halo. Gao et al. (2011) confirm that the fractional mass in subhaloes, 
)J>N(> /x), gently increases until [i « 0.01 and falls steeply above. The peak in the 
fractional mass increases from w 0.1% to « 1% from galaxy- to cluster-sized host 
haloes. 

1.4 Dark-matter annihilation? 

Recent observations and surveys of high-energy gamma rays with imaging Chercnkov 
telescopes and with the Large Area Telescope (LAT) on-board the Fermi satellite have 
resulted in a variety of theoretical investigations of galaxy clusters and interesting 
inferences. So far, galaxy clusters have escaped detection both by Cherenkov imaging 
and by Fermi-LAT. Both types of telescope cover very wide ranges of gamma energies. 
While the High-Energy Stereoscopic System (H.E.S.S.) is sensitive to gamma photons 
with energies between ~ 10 GeV to ~ 10 TeV, Fermi-LAT shows the highest response 
at smaller gamma energies between ~ 1 GeV and ~ 100 GeV. 




(16) 
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After one year of operation, the non-detection with Fermi-LAT of a gamma signal 
from nearby galaxy clusters led Yuan et al. (2010) to conclude that cluster substruc- 
tures should have masses of at least M > 10 2 . . . 10 3 M©. The argument behind this 
conclusion is straightforward: Substructures in clusters are expected to be substan- 
tially denser than the surrounding, smoothly distributed cluster matter, and are thus 
expected to emit a good fraction, if not most of a cluster's hypothetical gamma signal. 
The non-detections of gamma rays from nearby clusters then imply an upper limit to 
the lumpiness of cluster matter. Assuming a mass function for the cluster substruc- 
tures, this upper limit translates to a lower limit of substructure masses, since the mass 
function should increase quite steeply towards decreasing mass. Sanchez-Conde et al. 
(2011) used essentially the same line of reasoning to argue that nearby galaxy clusters 
may be more efficient probes of dark-matter annihilation than local dwarf galaxies 
when viewed with imaging Cherenkov telescopes. Conversely, Pinzke et al. (2011) ar- 
gued that, should clusters contain substructures down to the Earth's mass, models 
of dark-matter annihilation through leptonic channels (so-called leptophilic models) 
would be ruled out by the lack of detections. 

Similarly, the cross section for dark-matter annihilation was constrained from above 
by the gamma non-detections. Ando and Nagai (2012) used 2.8 years of Fermi-LAT 
data taken from 49 nearby clusters to derive (av) < 10 _ cm s _ at 10 3 GeV for 
the velocity-averaged annihilation cross section without taking into account that bary- 
onic energy dissipation may contract the dark-matter distribution, and argued that 
baryonic contraction could lower this upper limit by about an order of magnitude. 
From three years of Fermi-LAT data and gamma non-detection from eight nearby 
galaxy clusters, Huang et al. (2012) derived a lower limit on the dark-matter decay 
time of r > 10 26 s and an upper limit on the dark-matter annihilation cross section of 
(av) < 10 _22 cm 3 s _1 at gamma energies near 10 3 GeV. These constraints tighten by 
~ 3 orders of magnitude once cluster substructures are taken into account. Abramowski 
et al. (2012) placed an upper limit on the dark-matter annihilation cross section of 
(av) < 10 _23 cm 3 s _1 at 10 3 GeV, possibly lowering by > 2 orders of magnitude due 
to substructure, from the non-detection of the nearby Fornax cluster in 14.5 hours of 
observation with H.E.S.S. 

These conclusions are well consistent with each other and show that upper limits 
on the annihilation cross section as tight as (av) < 10~ 25 cm 3 s _1 can be inferred at 
gamma energies near ~ 1 TeV from the gamma non-detections if galaxy clusters are 
substructured as expected from simulations. Interestingly, Gao et al. (2012) deduce 
from highly resolved simulations not only that the dark-matter annihilation signal 
from galaxy clusters should be stronger than from galaxies because of the substruc- 
tures, but also that it should emerge preferentially from large cluster-centric radii since 
substructures are likely to be tidally destroyed near cluster cores. They identify the 
Coma and Fornax clusters as most promising candidates for gamma detections in future 
observations. 

These results should be read with some caution, however. Firm predictions of the 
gamma emission expected from galaxy clusters require that the abundance and the 
central density of substructures with very low masses be accurately known, possibly 
down to stellar- or even lower-mass structures. This is well outside the reach of even the 
best-resolved numerical simulations today. Current inferences from the non-detection 
of gamma rays from clusters are thus based on substantial extrapolations of simulation 
results towards masses well below their resolution limit. 
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Summarising the discussion so far, we see that galaxy clusters, which are at the 
massive end of the spectrum of gravitationally bound cosmic structures, formed latest 
during the cosmic evolution if structure formation proceeded in a bottom-up way, 
as expected if the dark matter is cold. In such a scenario, the most massive objects 
must always be a dynamically active population since they are the most recent to have 
formed and did at least in part not have the time to relax. Small, lower-mass structures 
formed earlier, which means that the more massive objects must at least partially have 
been assembled by merging smaller units. Even though it is expected that part of the 
small structures that finally end up as parts of a galaxy cluster get disrupted, clusters 
should host sublumps with a broad continuum of masses. The mass fraction contained 
in the sublumps is expected to reach the per cent level, with subhaloes of approximately 
one per cent of the cluster mass contributing most. Tidal disruption of subhaloes near 
cluster cores renders the spatial subhalo distribution substantially flatter than that of 
the dark matter. The non-detection of gamma rays from clusters puts firm upper limits 
on the dark-matter annihilation cross section and on the decay time of the dark-matter 
particles. 



Fig. 4 Haloes made of cold dark matter are expected to contain rich substructure on all scales. 
(Image extracted from the Millennium simulation, Springel et al. 2005) 



1.5 Attempts at understanding universality 

As we have described in some detail above, clusters as well as their substructures 
are expected to have universal density profiles according to numerical simulations. 
Despite many attempts to understand which physical processes shape this profile in 
absence of an equilibrium state, it is probably fair to say that a fundamental reason 
has not been found yet. Two main theoretical concepts have so far been discussed 
as possible explanations. One relies on the assembly of haloes from smaller haloes 
and the importance of the merging processes and rates involved (Syer and White 
1998; Subramanian et al. 2000). If n is the slope of the dark-matter power spectrum 
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near the scale corresponding to the halo mass considered, then halo density profiles 
should be power laws with slope —3(3 + n)/(5 + n). For n w —2 as appropriate for 
cluster-sized haloes, slopes near —1 can be achieved, as simulations find them near the 
cores of dark-matter haloes. Since the universal profile is found for haloes with masses 
spanning several orders of magnitude, haloes of galaxy mass for which n approaches —3 
should have natter cores than haloes of cluster mass, which is contrary to other claims 
based on simulations (Jing and Suto 2000). Moreover, halo mergers are important for 
this scenario, which are frequent in any bottom-up hierarchical scenario of structure 
formation such as cold dark matter (see Fig. 5). Simulations undertaken with hot dark 
matter, in which structure formation proceeds from the top down and mergers are rare, 
however, produce haloes with the same universal density profile, but with an inverted 
concentration-mass relation. Mergers can therefore not be essential for the formation 
of the universal profile (Moore et al. 1999; Wang and White 2009). 
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Fig. 5 Amount of mass accreted by a halo of mass 5 X 10 14 Mq per unit logarithmic scale- 
factor as a function of the ratio between the masses of the main halo and the merging haloes. 
At redshift zero, most of the mass is accreted in lumps having about half of the main-halo 
mass. At higher redshift, most of the mass comes in smaller pieces. 



Another extensively discussed concept relies on the secondary-infall model (see 
Ascasibar et al. 2007, for a recent example), which essentially rests on the fact that 
the so-called radial action is an adiabatic invariant. Generally, for a mechanical system 
with a Hamiltonian function varying slowly compared to the orbital time, the quantity 

P ' dq (17) 

is invariant in the adiabatic approximation, where q and p are the vectors of the 
generalised coordinates and their conjugate momenta. Particles of mass m on circular 
orbits in the potential of a spherically-symmetric mass M(r), for example, have circular 
velocities 

, / GM(r) 
Mr) = y ^ ■ (18) 

For their orbits, therefore, the quantity 

I = 2irm^GM(r)r (19) 
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is invariant while the potential changes adiabatically. For matter orbiting in exist- 
ing potential wells that deepen slowly by accreting matter from their environment, 
the square root of the radius times the mass is adiabatically invariant. While this as- 
sumption yields density profiles approximating the universal profile, the scatter due to 
varying initial conditions and angular momenta is large. Besides, orbital times are of 
order 

3/2 

t ^ s 4 w ' (20) 

while the time scale for changes in the potential may be as fast as the free-fall time 
scale, 

T « ' (21) 

which is of the same order as the orbital time. This sheds doubt on the adiabatic 
assumption at least in such circumstances when rapid mergers change the potential 
on a short time scale comparable to the free-fall time scale. Other discussions begin 
by postulating a certain phase-space structure of the dark-matter halo particles, for 
example a power law in the ratio p/a 3 between the density and the velocity dispersion 
a (Austin et al. 2005). Among others, the Jeans equation then admits constrained 
solutions which may come near the universal profile, but then the origin of the initial 
phase-space structure remains unexplained. 

Of course, this brief discussion cannot do justice to the many efforts to explain the 
fundamental problem why a physical process without an equilibrium state should reveal 
universal behaviour, but it may serve to illustrate why the mere existence of a universal 
profile, quite irrespective of its details, raises conceptually important questions that 
have so far escaped a fundamental explanation. 



1.6 Halo concentrations and triaxiality 

Comparably fundamental may be the concentration parameter which controls the shape 
of the universal profile. Simulations show that it should decrease with mass quite slowly, 
c oc M -0 ' 1 or similar, and more steeply on redshift. ft should be mentioned that there 
is at least one simulation (Prada et al. 2012) claiming that the concentration-mass 
relation be non-monotonic, increasing towards high masses. At the time of writing of 
this review, it is yet unclear why this result disagrees with virtually all other simulations 
of dark-matter haloes, and whether it may be seen as an artefact of interpretation. 

The decrease of the concentration with mass is generally explained as a reflection 
of halo formation time. As mentioned above, it is assumed that halo cores preserve the 
mean cosmic density at the time of their formation, which at least qualitatively repro- 
duces the trend seen in simulations of hierarchical models of structure formation such 
as cold dark matter: Massive haloes form later and should thus have lower concentra- 
tions. This is supported by the finding that the concentration-mass relation is inverted 
in top-down formation scenarios occurring in warm dark-matter models (Wang and 
White 2009). 

As is detailed elsewhere (see the review by Ettori et al. in this volume), some mea- 
sured concentrations disagree strongly with those expected from the concentrations- 
mass relations found in numerical simulations (Broadhurst et al. 2005, 2008; Comerford 
and Natarajan 2007; Hennawi et al. 2007; Lemze et al. 2009; Umetsu and Broadhurst 
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Fig. 6 Left panel: Examples of concentration-mass relations at redshift zero in the mass range 
of galaxy cluster according to different algorithms relating the formation time to the halo core 
density. Right panel: Evolution of the concentration of haloes with mass 5.0 X 10 14 M Q with 
redshift according to the same algorithms. (Bullock et al. 2001; Eke et al. 2001; Navarro et al. 
1997, 1996; Seljak 2000; see Dolag et al. 2004; Duffy et al. 2008; Gao et al. 2008; Maccio et al. 
2008, 2007; Neto et al. 2007; Shaw et al. 2006; Zhao et al. 2009 for some other examples) 



2008; Umetsu et al. 2010). In particular, concentrations derived from strong gravita- 
tional lensing tend to be substantially and significantly higher for massive clusters, 
sometimes more than twice as high. For massive clusters, concentrations around c = 5 
should be found, as Fig. 6 shows. At fixed mass, concentrations of simulated haloes 
have a log-normal distribution with a standard deviation of w 0.2, corresponding to 
typical relative concentration fluctuations of « 20 per cent around the mean. Concen- 
trations as high as those routinely found in strongly-lensing clusters are difficult to 
understand since they can hardly be accommodated with the standard /tCDM cos- 
mological model. Baryonic physics and feedback processes in cluster cores may affect 
the expected concentrations, though, as discussed by Mead et al. (2010). They find 
that higher concentrations occur if baryonic cooling is taken into account, while AGN 
feedback brings the concentrations back near the range expected without baryons. 

However, a substantial concentration bias can be expected in galaxy-cluster sam- 
ples selected for their strong-lensing effects. Selecting strongly gravitationally lensing 
cluster-sized haloes in simulations, Meneghetti et al. (2010a) find that such clusters 
have concentrations tend to have concentrations higher by m 50 % than average. Tak- 
ing this bias into account, Umetsu et al. (2011) find that their high concentrations 
derived from strong and weak lensing in four similar clusters maybe compatible with 
theoretical expectations. 

However, high concentrations in clusters are not an unanimous result. Some strong- 
lensing studies, sometimes even of the same clusters from which high concentrations had 
been inferred, find concentrations in the expected range (Halkola et al. 2006; Limousin 
et al. 2008). In a strong-lensing analysis of 12 clusters at redshift z > 0.5 from the 
Massive Clusters Survey (MACS), Sereno and Zitrin (2012) find good agreement of 
the measured concentrations with theoretical expectations if the possible triaxiality of 
the cluster mass distribution is taken into account, even though their effective Einstein 
radii are still substantially higher than expected for ylCDM. To further complicate 
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the picture, concentrations inferred from weak gravitational lensing generally seem 
to reproduce the expected concentrations well. To give just two recent and possibly 
extreme examples, Lerchster et al. (2011) find a concentration of c = 4.0^.!, 4 from a 
weak-lensing analysis of the cluster XMMU J1230. 3+1339 at z = 0.975. Okabe et al. 
(2010) analyse the concentration-mass relation with weak gravitational lensing in 30 
X-ray luminous clusters. They find a fairly steep dependence of the concentrations on 
mass, c oc M _0 - 40±0 - 19 j DU t the inferred concentration of c = 3.48^}'®5 for clusters 
with a virial mass near M = 10 15 /i _1 M© are just as expected. Yet another result 
is found by Mandelbaum et al. (2008), who find concentrations lower than expected 
in a large sample of haloes from a wide mass range analysed by weak lensing. Using 
numerical simulations, Bahe et al. (2012) find that the concentrations expected to be 
measured from weak lensing should be slightly biased low. 

In contrast to gravitational lensing, X-ray observations (e.g. Buote et al. 2007; 
Ettori et al. 2011, 2010; Gastaldello et al. 2007; Pointecouteau et al. 2005; Schmidt 
and Allen 2007) and interpretations of galaxy kinematics in clusters (e.g. Rines and 
Diaferio 2006; Wojtak and Lokas 2010) typically find good agreement between mea- 
sured and expected concentrations, in particular for clusters considered relaxed, albeit 
with considerable scatter. Detailed analyses of cluster cores based on X-ray data or 
on lensing data combined with the stellar dynamics in the central cluster galaxy find 
a large spread of density profiles, ranging from noticeably steeper than NFW (Lewis 
et al. 2003; Buote and Lewis 2004) to considerably flatter (Sand et al. 2004; Newman 
et al. 2013; Ettori et al. 2002). Clearly, no final word can yet be said on the central 
density profiles of clusters. An outstanding example of a cluster showing an NFW den- 
sity profile is Abell 2261 (Coe et al. 2012), observed within the CLASH programme 
(Postman et al. 2012). 

Many recipes exist giving empirical descriptions of the concentration-mass rela- 
tion. A selection of them is displayed in Fig. 6. Typically, they define a halo-formation 
redshift as the redshift when a certain fraction of the final halo mass was assembled, 
and relate the halo core density to the mean cosmic density at that redshift. Several 
of these recipes (Navarro et al. 1997; Bullock et al. 2001; Eke et al. 2001; Neto et al. 
2007; Seljak 2000; Zhao et al. 2009) are found to fit the mass-concentration relations 
of simulated haloes quite well, at least in certain wide mass and redshift ranges. How- 
ever, the formation redshifts defined that way have to be very high for the simulated 
concentration-mass relation to be explained. Apparently, not more than a few per cent 
of the final halo mass need to be assembled when the final core density of the forming 
halo is imprinted. At the very least, this is an astonishing finding that should ultimately 
be understood. 

Another difficulty related to the halo concentrations has to do with the so-called 
halo model of non-linear cosmic structure formation. It assumes that all cosmic struc- 
tures are composed of haloes whose two-point correlation defines the power spectrum on 
large scales, while the small-scale power spectrum is given by the mass auto-correlation 
inside individual haloes. This being so, the internal density structure of individual 
haloes is an important ingredient of the halo model. However, for reproducing the 
non-linear, small-scale power spectrum found in numerical simulations, concentration- 
mass relations typically have to be assumed which are substantially steeper functions 
of mass than found by direct measurement in the very same simulation data. It seems 
that this substantial discrepancy can be solved by taking into account that haloes are 
substructured and thus contain a hierarchy of lower-mass haloes (Giocoli et al. 2010a). 
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It can be shown analytically (Doroshkevich 1970) that cosmic structure formation 
by gravitational collapse from an initially Gaussian random field has to proceed in an 
anisotropic way. The probability for any pair of eigenvalues of the so-called Zel'dovich 
deformation tensor to be equal vanishes identically, hence structures have to evolve to- 
wards triaxial bodies at least well into the mildly non-linear regime. It is thus natural 
to expect that the shape of dark-matter haloes should approximate triaxial ellipsoids 
characterised by two axis ratios, called the ellipticity and the prolaticity, for the sta- 
tistical distribution of which any structure-formation scenario based on an Gaussian 
initial conditions makes specific predictions. Modelling of strong gravitational lensing 
(see there) shows that triaxial haloes can better fit the gravitational-lensing effects 
simultaneously with the X-ray emission of galaxy clusters, and the axis ratios found 
that way are generally in a reasonable range. 

The brief review given so far shows that the core structure of galaxy clusters poses a 
number of substantial problems, some of which are likely to reach far into fundamental 
physics. Cast into questions, these problems can be enumerated as: (1) How can self- 
gravitating systems like dark-matter haloes form a universal density profile despite the 
absence of an equilibrium state? (2) Do the density profiles of galaxy clusters adopt 
the shape found in numerical simulations? (3) How concentrated is the dark matter 
in their cores? (4) Are galaxy clusters as substructured as expected for haloes formed 
from cold dark matter? (5) Are clusters as triaxial as they are expected to be? 

Questions (1), (2) and possibly (3) address the fundamentally unsolved problem 
of what controls the non-equilibrium dynamics in of self-gravitating systems. Question 
(3) may be, question (4) certainly is related to the properties of the yet hypothetical, 
presumably weakly interacting, possibly annihilating dark-matter particle. Question (5) 
pertains our understanding of cosmic structure formation as being due to gravitational 
collapse from a Gaussian random field of density fluctuations. 



2 Baryonic physics in cluster cores 

In the centres of galaxy clusters, the mass fraction of baryons becomes significant. 
Whether it is large enough to substantially alter the properties of the dark matter 
"backbone" we do not know, but it is clear that one has to consider the influence of 
processes such as radiative cooling, heating and other feedback processes from baryons 
to dark matter when simulating cluster physics. Various processes act on scales that 
are orders of magnitude apart, which renders the simulation of the physics in cluster 
cores complicated. This is the reason why currently only approximate simulations of 
the conditions in cluster centres are available. 

Similarly, the interpretation of observations is not easy. An outstanding example 
of the complications encountered in cluster centres is the Perseus cluster (see Fabian 
et al. 2011, and references therein). This cluster can be studied in great detail due to its 
proximity and could thus be considered a prototype for the overall population. In this 
cluster, which on large scales is a regularly shaped, relaxed object, bubbles, ripples and 
shock fronts abound in the cluster core. There is a rich variety of phenomena including 
a complex gas temperature and metallicity structure, extended radio emission coming 
from the X-ray bubbles to filamentary structures that are bright in H-a, but also visible 
in the X-ray regime. In this system, it is also apparent that projection effects have to 
be taken into account when trying to understand the complicated interactions. 
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It has been shown (e.g., Voit et al. 2003) that the quantity 




has the properties of an entropy measure ("astronomer's entropy") in the sense that it 
is constant during adiabatic processes. Here T is the temperature and n e is the electron 
density, which is proportional to the gas density. When this quantity is plotted as a 
function of cluster radius and compared to results from non-radiative simulations (e.g., 
Pratt et al. 2010), the deviations are strongest in the cluster core in the sense that the 
observed K profiles lie above the numerical expectations. This is attributed to the effect 
of non-gravitational processes, or to processes that have acted before ("preheating") 
or during the formation of the structures. It serves to show the influence of non- 
gravitational processes on cluster-sized structures. Already Ponman et al. (2003; see 
also the review by Giodini et al. in this volume) found that when computing the 
measure K inside ten percent of the virial radius for several dozens of systems, their 
K-T relation did not match the simple power-law K oc T expected for systems which 
share the same average gas density, for example if they formed at similar redshifts. 
They concluded that non-gravitational processes are in fact at work not only for the 
coolest of their studied systems, but also for systems with temperatures of 5 keV and 
above. 

There is evidence for feedback processes in the centres of the most massive galaxy 
clusters. Here the gas densities of the order of ~ 1 particles cm -3 theoretically imply 
large cooling rates of several hundreds of solar masses per year in some extreme cases 
(e.g. Fabian 1994; Allen 2000). Although steep profiles of the temperature as a function 
of radius (projected and de-projected) have been observed in many cases (e.g. Allen 
et al. 2001; De Grandi and Molendi 2002; Vikhlinin et al. 2005; Pratt et al. 2007) that 
lead from 5 to 10 keV at r25oo for the most massive clusters down to around 1 keV in 
the cores of the clusters, no evidence has been found from radiative cooling that large 
quantities of gas are cooling below this value (e.g. Peterson et al. 2003). X-ray cavities 
reminiscent of those observed in the Perseus cluster have been found in many of these 
systems (e.g. Allen et al. 2006; Blanton et al. 2011). This suggests that heating through 
AGN feedback is an important mechanism in cluster cores. The exact way by which 
these different processes acting on very different size scales are balancing each other 
is still being studied (e.g. Vazza et al. 2012). Also, other processes such as thermal 
conduction (e.g. Voigt et al. 2002) or mixing have been studied (Fabian et al. 2001) 
and are still under scrutiny. 

A further process that has been repeatedly invoked in the interpretation of data 
is adiabatic contraction. The process leads to the contraction of the dark cluster halo 
due to the cooling and condensation of gas in cluster centres (e.g. Gnedin et al. 2004). 
Although such models can be fitted to data, it is not clear yet whether this process is 
at work for galaxy clusters. Despite all the described "baryonic impediments", there 
is an intermediate region with radii of the order of the NFW scale radius where the 
effects of the central dark matter cusp are measurable, but where the feedback effects 
are not (yet) affecting the measurement significantly (e.g. Pointecouteau et al. 2005; 
Schmidt and Allen 2007; Vikhlinin et al. 2006). Here, X-ray emission from the bary- 
onic component together with the assumption of hydrostatic equilibrium can be used 
to constrain the mass profile or the mass-concentration relation of galaxy clusters. The 
X-ray results on the cores of galaxies overall suggest a good agreement with the results 
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from numerical simulations in terms of the mass profile (e.g., Allen et al. 2001; Arabad- 
jis et al. 2002; Lewis et al. 2003; Voigt and Fabian 2006; Schmidt and Allen 2007; Host 
and Hansen 2011, and references therein). The determination of the amplitude of the 
mass-concentration relation is still somewhat uncertain and very likely affected by the 
large intrinsic scatter that is also observed in the simulations (Buote et al. 2007; Ettori 
et al. 2010, 2011; Schmidt and Allen 2007). 

3 Combining cluster observables 

How can a consistent picture of galaxy clusters be constructed that incorporates all 
kinds of information these objects provide? Of course, it is a viable approach to re- 
cover e.g. a cluster's radial mass profile by fitting a parametric model to X-ray and 
lensing data, as innumerable studies have shown. As the number of available empirical 
constraints on cluster mass distributions is increasing due to impressively improving 
observational data obtained in a multitude of channels, the question may be raised how 
all cluster observables may be combined on a common basis. 

Current cluster observables are their gravitational lensing effects, their X-ray emis- 
sion, the kinematics of their member galaxies, and the thermal Sunyaev-Zel'dovich 
effect. How can these different observables be combined on a common ground? The 
final section of this review sketches a possible approach based not on the mass or the 
mass density, but on the gravitational potential. 

3.1 Strong lensing and weak shear 

In the approximation relevant for our purposes, gravitational lensing is not affected 
by the physical state of the matter and energy distribution inside the lensing object. 
It thus provides perhaps the most reliable techniques for analysing the actual mass 
distribution in clusters (see e.g. Bartelmann 2010, for a review). Problems are generated 
by ambiguities due to a priori unknown degrees of asymmetry, line-of-sight projection 
effects and the limited range of constraints in particular by strong gravitational lensing. 

Strong lensing is being exposed in detail elsewhere in this volume. The main state- 
ments explained there are as follows: Strong lensing occurs with sources close to caustic 
curves, their images occur near critical curves. Two types of critical curve can be dis- 
tinguished by the direction into which images are preferentially stretched. Tangential 
critical curves constrain the total projected amount of mass enclosed by them, albeit 
that mass is scaled by the critical surface mass density which depends on cosmology and 
sometimes unknown redshifts. Radial critical curves, which are typically surrounded 
by the tangential critical curves, constrain the local slope of the density profile. 

The most prominent phenomena of strong gravitational lensing are tangential and 
radial arcs. The latter are difficult to find since they occur near cluster cores, where 
they may be overshone by the central cluster galaxies. Arcs trace or straddle the critical 
curves they are caused by. Multiple images, which are much more frequent than arcs, 
also constrain the locations of the critical curves. 

One might think that strong lensing should constrain the core structure of galaxy 
clusters quite well. However, this is not the case if model-independent constraints are 
sought. Since all phenomena of strong lensing rely in some way or another on the 
occurrence and the locations of critical curves, all the information they provide concern 
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the local curvature of the effective lensing potential. Strong-lensing observables can 
fix parametric mass models quite tightly, but they rarely allow model-independent 
statements on their own. 

Weak gravitational lensing (cf. Bartelmann and Schneider 200f ) weakly distorts the 
images of background galaxies. To lowest order, weak lensing imprints an ellipticity on 
background sources characterised by the local strength of the gravitational tidal field or 
shear 7 and the scaled surface-mass density or convergence k. Both k and 7 are linear 
combinations of the effective gravitational potential ip. The shear has two components, 

71 = \ (dfip - dli^ , 72 = did 2 ip , (22) 

which can conveniently be combined to form the real and imaginary parts of the com- 
plex shear 7 = 71 + i72- The convergence k is half the Laplacian of ip, 

*=\ (dU + 6$1>) • (23) 

To lowest-order approximation of weak lensing, idealised, circular sources appear as 
elliptical images. Let a and b be the semi-major and semi-minor axes of such an ellipse, 
then its ellipticity 

(24) 

a + b y ' 

is approximately given by the so-called reduced shear g rather than the shear itself, 

e »° := ih- (25) 

Circular sources are rare if they exist at all, but can be approximately created by 
averaging over sufficiently many background sources if their intrinsic shapes are random 
and uncorrelated. Thus, weak-lensing effects have to be inferred in a statistical way 
by measuring the net ellipticity after averaging over images. This immediately implies 
that weak lensing has an intrinsic resolution limit set by the mean number density n 
per solid angle of background sources. If N galaxies should be averaged over, they need 
to be collected within a finite solid angle whose radius 9 defines the angular resolution. 
The condition N > niv8 2 requires 

/ AT 

(26) 

Typical numbers are n w (30 . . . 40) arcmin" 2 and N « 10, giving 6 > (17... 20)". 
The resolution limit of weak-lensing observations thus turns out to be of the same 
order as the typical angular extent of the critical curves. Strong-lensing observables 
thus constrain effective cluster potentials on scales at or below the angular resolution 
provided by weak-lensing observables. If they are to be combined, and they need to be 
in order to constrain the core structure of galaxy clusters in a model-independent way, 
a multi-scale approach is required. 

Such multi-scale approaches have been developed and applied (Bartelmann et al. 
1996; Bradac et al. 2006, 2005; Cacciato et al. 2006; Meneghetti et al. 2010b; Merten 
et al. 2009, 2011; Umetsu et al. 2012). In essence, they map the effective lensing poten- 
tial on a grid with a resolution adapted to the observed constraints, derive its lensing 
properties by discrete differentiation and vary the gridded potential values such that 
a suitable total likelihood is minimised. This likelihood contains terms from strong 
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and weak lensing. It expresses the expectation that the reconstructed lensing potential 
should recover the measured mean ellipticities and the location of the critical curves 
traced by any strong-lensing phenomena. Simulations show that such parameter-free, 
adaptive-grid methods recover the density profiles and the total cluster masses very 
well (Meneghetti et al. 2010b). 

Multiple images constrain the deflection angle near critical curves, and thus the 
gradient of the lensing potential. As described elsewhere, critical curves are located 
where the determinant of the Jacobian matrix vanishes, which essentially provides a 
constraint on the matrix of second derivatives (the Hessian matrix) of the lensing po- 
tential. The gravitational shear measured from elliptical distortions constrains another 
combination of second-order derivatives of 4>- The different degrees of derivatives enter- 
ing these constraints are conceptually important since the trace the lensing potential on 
different scales: The higher the derivative of the potential is that underlies an observ- 
able, the smaller are the scales that the respective observable is sensitive to. Aiming at 
the core structure and the amount of substructures in clusters, combining observables 
which are sensitive on different scales is conceptually important. 

3.2 Gravitational flexion 

Can third-order derivatives of the lensing potential be used? Their effect can straight- 
forwardly be derived in the following way (Goldberg and Bacon 2005; Bacon et al. 
2006). Combine the derivatives with respect to the two spatial directions into one 
complex- valued differentiation operator d and its complex conjugate d* , 

d:=d 1 + id 2 , d* := 9i - id 2 ■ (27) 

The convergence and the shear can then simply be written as 

K = \ertoi>, 7 = ^ 2 v>. (28) 

Further differentiation defines the two flexion components F and G, 

F = dn , G = d~f , (29) 
with both F and G being complex numbers. Expressed in components, we have 

F\ = UF = ^ (Vm + V122) , F 2 = <3F=^ (V112 + V-222) (30) 

and 

Gi = KG = i (V111 - 3^122) , G 2 = 3G = ^ (3^112 - ^222) • (31) 

From the construction of the flexion components, it is quite easily seen that F and 
G create images with one- and three-fold rotational symmetry from circular sources, 
respectively. Such images are said to have spin-1 and spin-3. While the spin-1 image 
distortion due to the F flexion causes a centroid shift which may be hard to recognise 
observationally, the spin-3 distortion due to the G flexion imprints a characteristic 
triangular shape. 

Faint galaxies tend to appear intrinsically elliptical. The gravitational shear, adding 
elliptical distortion, thus has to be detected on top of the intrinsic galaxy ellipticity. In 
contrast, since there should be little or no intrinsic triangular morphology in galaxies, 
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Fig. 7 Illustration of gravitational lonsing effects on a circular source. Left column: shear; 
centre column: F flexion; right column: G flexion 



the triangular distortion due to the G flexion should be measurable with much lower 
intrinsic noise than the shear, and much less systematic uncertainty than the centroid 
shift caused by the F flexion. The signal-to-noise ratio specifically of the G flexion may 
thus be comparable to that of the shear, even though the signal itself is expected to 
be weaker. 

To give an order of magnitude, consider the simple lens model of an isothermal 
sphere, whose intrinsic, three-dimensional density profile is proportional to r~ 2 . Its 
convergence and the absolute value of the shear are given by 



(32) 



where 6*e is the Einstein radius. Clearly, k and 7 are dimension-less. The F and G 
flexions have the amplitude 



\m\ = ^ = \G(9)\ 



(33) 



and the dimension of an inverse angle. Their signal, measured across an image of 
characteristic angular scale 88, is thus estimated to be 



\F(0)\se = e ^j = \G(e)\59, 



(34) 



lower by a factor 58/8 than the shear signal. This factor is should fall between 0.1 and 
0.01 for typical weakly-lensed images. 

It is not clear yet whether it will be possible to reliably measure flexion in galaxy 
clusters. Some pioneering attempts have been made which found hints at F rather 
than G flexion (e.g. Okura et al. 2007; Goldberg and Leonard 2007; Leonard et al. 
2007; Okura et al. 2008; Schneider and Er 2008; Pires and Amara 2010; Velander et al. 
2011; Leonard et al. 2011; Cain et al. 2011). However, flexion can only be measured if 
a reliable estimate for the shear already exists because any spin-1 or spin-3 quantity 
that can be derived from observations of galaxies not only picks up the flexion signal, 
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but combinations of flexion and shear (Viola et al. 2012). Considerable work on shape 
measurements of galaxies and their interpretation still needs to be done before flexion 
can be considered measurable. If flexion can be determined, however, it promises to 
add small-scale information on the internal structure of galaxy clusters that the shear 
alone cannot resolve. 

The essential message from this subsection is two-fold: Gravitational lensing primar- 
ily provides information on the second-order derivatives of the gravitational potential. 
This is necessarily so because of the equivalence principle: It is not the gravitational 
field that can be invariantly measured, but the gravitational tidal field, which is in the 
weak-field limit represented by the Hessian of the Newtonian potential. Should flexion 
become reliably measurable, it would allow quantifying changes in the tidal field, sup- 
plying smaller-scale information than the gravitational shear, the magnification and 
critical curves can provide. Lensing observables thus constrain the projected gravita- 
tional potential by its curvature and possibly higher-order derivatives. Gravitational 
lensing returns an estimate for the projected Newtonian potential up to harmonic 
functions. 



3.3 X-ray emission 



X-ray emission by thermal bremsstrahlung can also be analysed in terms of the grav- 
itational potential. To illustrate this statement, assume that the X-ray emitting gas 
with density p, pressure P and temperature T is in hydrostatic equilibrium with the 
gravitational potential <P, satisfying 



VP 



= -V<2> . 



(35) 



Assuming a polytropic stratification of an ideal gas with an effective polytropic index 

7, 

P = Po(^)\ (36) 

the hydrostatic equation (35) can easily be solved to give the gas density in terms of 
the gravitational potential, 



Po 



7 - 1 P0 
7 Pq 



1/(7-1) 



(37) 



where the potential <£rj enters as an integration constant that can be set to ensure 
p = at a finite radius from the cluster core. By the ideal gas equation, (37) implies 
the relation 

7- 1 PO 



T 



7 Po 



(* - $o) 



(38) 



The fiducial pressure, density, and temperature, Po, po and Tq, can be fixed at an arbi- 
trary cluster radius, e.g. at the centre. Equations (37) and (38) imply the (bolometric) 
bremsstrahlung emissivity 



Jx = Jx.o 



7 - 1 PO 
7 Po 



(39) 
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with an exponent 



+ 1 = ^- (40) 



7-1 2 2(7-1) 



For 7 « 1.2, which may be appropriate for an estimate, a « 10.5: The X-ray emissivity 
depends very steeply on the gravitational potential. The exact value of the constant 
jx o is known but irrelevant for our purposes. 

The relation (39) that we have derived between the X-ray emissivity by thermal 
bremsstrahlung and the Newtonian potential suggests a path from the X-ray surface 
brightness Sx to the projected gravitational potential, thus connecting X-ray emis- 
sion and gravitational lensing on a common ground. By means of Richardson-Lucy 
deprojection, an iterative scheme can be constructed allowing to recover the three- 
dimensional X-ray emissivity jx from Sx- The Richardson-Lucy deprojection rests on 
Bayes' theorem and returns, usually with a regularisation term suppressing noise over- 
fitting, a three-dimensional function from its two-dimensional projection, provided a 
known projection kernel. 

Notice that we have not applied any symmetry assumptions on the X-ray emissiting 
cluster so far. The integration of the hydrostatic equation leading to (37) and thus also 
to (38) is possible without adopting a particular symmetry. This has to be given up 
now before Richardson-Lucy deprojection can be applied because the projection kernel 
depends on the symmetry of the three-dimensional shape assumed for the body of the 
hot intracluster gas. Assuming a reasonable symmetry (which by no means has to be 
spherical), Richardson- Lucy deprojection returns an estimate for jx from the observed 
S x . 

Normalising the result, and raising it to the power of 1/a, (39) allows us to derive 
the scaled potential 

£> ( <2>_<2> o) = ^(22L.) 1/a . (41) 

Po V 0> 7-1 \jx,oJ ( ' 

The exponent a needs to be known, but can be derived via the polytropic index 7 
which can be determined from the temperature profile. Notice that the expression on 
the left-hand side of (41) is dimension-less since the potential is scaled by the squared 
thermal velocity Po/po, which can again be estimated from the measured temperature 
and an assumed mean gas-particle mass. The scaled potential can then be projected 
along the line-of-sight to obtain an independent estimate of the projected Newtonian 
potential, which can directly be combined with the estimate of the projected potential ip 
provided by lensing. First tests of this algorithm assuming spherical symmetry returned 
very promising results. 



3.4 Thermal Sunyaev-Zel'dovich effect 

The same procedure developed above can now also be applied to observations of the 
thermal Sunyaev-Zel'dovich effect. This effect is quantified by the Compton-j/ param- 
eter 

V = ^vrne , (42) 

where the relevant quantities are the electron number density n c and the gas temper- 
ature T and the other quantities are constants. Since the electron number density is 
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determined by the gas density, we can use the results (37) and (38) obtained above to 
write 



(43) 



v = yo 

with the exponent 



)9 = — Kr + 1 = w 8.3 . (44) 

7 — 1 7—1 

Now, the same algorithm can be applied as for the X-ray surface brightness: The 
observed, line-of-sight integrated Compton-y parameter Y can be Richardson-Lucy 
deprojected to yield the three-dimensional Compton-y parameter, which can the be 
converted to the (thermally scaled) gravitational potential and projected, 

£(*-*0) = -^ L 1 ( JL ) 1/t> ■ (45) 
Pa 7 - 1 V VO J 

Again, first tests show that this procedure returns reliable results. 

The main result of the last two subsections is that the observables provided by the 
hot intracluster gas can directly be converted to constrain the projected Newtonian 
potential. This allows combining all related observables - strong lensing, weak lensing 
possibly including flexion, X-ray emission and the thermal Sunyaev-Zel'dovich effect - 
into one x 2 function returning a joint estimate for the projected cluster potential. 



3.5 Galaxy kinematics 

One observable remains, which is related to the kinematics of the cluster galaxies. 
Assuming spherical symmetry in the Jeans equation, we can write 

d(no 2 r ) 2p{r). 2 . d<P 

__ + _ (nar) = _ n _, (46 ) 
s the spatial numbt 

sion, and 



where n is the spatial number density of the galaxies, a 1 is their radial velocity disper- 



P{r) = l--f (47) 

is the usual anisotropy parameter of the velocity ellipsoid. We shall later assume for 
simplicity that the number density n is proportional to the matter density p. 

Instead of the usual approach, which would consist in solving (46) to arrive at a 
fairly involved integral relation between (noy) and d r $ for given /3(r), we can now 
follow a different route, relating the effective "galaxy pressure" P ga j = no 2 - directly to 
the matter density p assuming a polytropic stratification with an effective polytropic 
index k, 

Pgal = P S ,lfi (£) =^gal,o(^) (48) 

Numerical experiments show that for a wide range of assumed matter-density profiles 
and anisotropy parameters, such effective polytropic indices can indeed be found. Then, 
(46) can be transformed into the equation 
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for the function 




with the exponent 



The equation (49) can now be solved in terms of q first, resulting in a Volterra integral 
equation of the 2nd kind relating q to the gravitational potential. Now, the line-of- 
sight integrated "galaxy pressure" is the observable line-of-sight velocity dispersion cry , 
weighted by the galaxy number density. Again, by means of Richardson-Lucy deprojec- 
tion, this observable can be converted to the three-dimensional "galaxy pressure" P ga l, 
thus providing the measureable source term in the Volterra integral equation for the 
gravitational potential <P. This can then be solved by a standard, iterative procedure. 
The three-dimensional potential obtained this way can finally also be projected along 
the line-of-sight to return the projected gravitational potential if>. 

The main result of this subsection is that also the observable galaxy kinematics 
can be turned into a constraint on the two-dimensional gravitational potential, given 
an assumption on the anisotropy parameter of the galaxy orbits. 

In conclusion, we have seen in the last section that all cluster observables can be 
combined on the grounds of the projected gravitational potential, jointly probing its 
structure on different scales with different resolution. First tests of these approaches 
have been undertaken with promising results. Details remain to be worked out, but it 
seems possible that this joint approach will return the tightest possible constraints on 
the gravitational potential, and thus on the density distribution, in clusters and their 
cores. 
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